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1 Introduction 


The computation of correlation functions is one of the major challenges in integrable 
systems. Multidimensional integral formulae have been derived for the infinite PHEl 
as well as the finite XXZ spin-chain at zero HE) and at finite temperature [B]. More 
recently it has been observed that some of these multidimensional integrals can be 
reduced to one-dimensional ones allowing for their explicit computation [3 IHl 1^ • In 
a series of papers unmiiiii] this reducibility was connected to a duality between the 
solutions of the quantum Knizhnik-Zamolodchikov (qKZ) equations of level 0 and level 
-4 which the correlation functions of the model ought to obey CSI. 

In this article we will refer to the subsequent work by Boos, Jimbo, Miwa, Smirnov 
and Takeyama on the correlation functions of the infinite XXZ spin-chain in, 


Hxxz = 9 X] 


j=-oo 


ctJuJ+i + + 


q + q 


-1 




q = e 


( 1 ) 


In there it is explained that the correlation functions are particular solutions of a reduced 
set of qKZ equations which can be formulated in terms of a special trace of a certain 
monodromy matrix. The latter formally resembles in its algebraic structure the mon- 
odromy matrix of the inhomogeneous six-vertex model on a finite lattice whose number 
of columns is determined by the number of operators in the correlation function. It is 
the special trace of this finite monodromy matrix which we will identify as an analytic 
continuation of the six-vertex fusion hierarchy in terms of Q’OP^i'citors. Before giving 
the technical details of the computation we state this result more explicitly. To this end 
we first briefly summarize the main outcome of the work m. 


1.1 Correlation functions and a generalized trace 

Consider the correlation functions of the elementary matrices {E^. )j acting on the 
site in the XXZ spin-chain and interpret 

p 

= (uac|(K_£^,gJi---(S_£p,gJp|uac)]J(-ej) (2) 

i=i 

with Sj,ej = ±1 as the vector components of a function hp = hp{Xi, ..., Ap) whose values 
lie in the tensor space V = Cu+ © Cu_. Here |uac) denotes the groundstate of 

the six-vertex model on an infinite lattice with inhomogeneity parameters {Aj}. In the 
homogeneous limit Xj 0 this becomes the groundstate of the Hamiltonian (^*. The 
vector valued functions hp are solutions of a set of reduced qKZ equations and according 
to m can be expressed as 

gf2p(Ai,...,Ap) P _ 

hp{Xi, ...,Xp) = - — - \_Y^e.,b 5 = v+®v-- e. = 2p+l - I (3) 

_^ 

*We follow here the conventions in M and consider in the massive regime, where the the groundstate 
is twofold-degenerate, only matrix elements between the same vectors. Notice further that unlike in the 
case of the XXX chain it is not known at the moment how to take the homogeneous directly in the 
ansatz of m 
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where Cl = Yli<j is a sum of operators of the form 

with Xij = Xi — Xj and uj,Lij being certain scalar functions. The operators W,W de¬ 
pending rationally on {g^*} are dehned through the aforementioned special trace of a 
monodromy matrix. For instance, setting {i,j) = (1,2) and Aq = (Ai + A2)/2 one 
considers the map l/® 2 (p- 2 ) y(g) 2 p 


2(Ai,..., Ap) — 


TrAi2 -^ 2(^02 — 1) • ■ ■ Lp{Xop — l)Lp(Aop) ■ ■ ■ T2(Ao2) 
[Al2]q 0^=3 


51,2^1,2 (5) 


which decomposes as 


p^p—2 (Ai, ..., Ap) — Ai 2 'pGp— 2(9 


Ai 


I +pGp-2{Q^^ , ■■■, 




( 6 ) 


The two terms appearing in the above sum determine the two operators on the right 
hand side in For further details and the precise relation between G, G and W, W 
we refer the reader to M- For our purposes it will be sufficient to focus only on the 
following object in Endl/®^^”^ 


t Tr^j^j — A 2 — 1) ■ ■ ■ Lp[X — Ap — l)Lp(A — Ap) • • • L2{X — A 2 ) . 
The operator Lj G Uq{sl 2 ) 0 End Vj- is a quantum group intertwiner and 

TrA = 


(7) 

( 8 ) 


denotes a linear function 


TVa,c : Uqish) 0 C[C, r '] - AC[C, C"'] 0 C[C, C"'] (9) 


which for A = n -|-1 yields the conventional trace over the quantum group representation 
7r0) of spin n/2, 

TrA=n-hi(a;) = Tr x, Vx G Uq{sl 2 ) ■ (10) 

7r(") 

Eor our definition of see (pn|) in the text below. In particular, one has for the special 
elements x = 1, g^ and the Casimir operator 


_ + q^- 

(g-g-i)2 


the following identities 


T1a,c 1 = A, TrA,cg™'‘ = 


qm _ g-m ^ 


Ti-a,c‘^ = 


C0C~^ 

(g-g-i)2 ■ 


( 11 ) 


( 12 ) 


We will now explain the main point of this article, namely, that there is an alternative 
expression for © which does not use the introduction of the abstract trace function 
® but relies on the two linearly-independent solutions to Baxter’s TQ equation for the 
inhomogeneous six-vertex model with twisted boundary conditions. 
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1.2 Baxter’s TQ equation and the six-vertex fusion hierarchy 

Consider the inhomogeneous six-vertex model on a finite lattice with length M = 2p—2 G 
2N, compare with ©• This model can be solved by finding solutions to Baxter’s TQ 
equation nnmainiiiH], 

M M 

t{\)Q{X) = Q{X + 1) n [A — Xm\q + Q{X — 1) n [A — Xm + l]q • (13) 

m=l m=l 

Here t denotes the six-vertex transfer matrix and Q is an auxiliary matrix. In terms of 
eigenvalues m is a difference equation of second order and hence will in general allow 
for two linearly independent solutions, say nainiiiii. However, not both of them 
always possess the same analyticity properties. For instance, for generic q and when M 
is even, as it is the case here, there is only one solution which can be expressed as a 
product of the form 

Q+(A) = n [A - = n ^^ Bethe roots, (14) 

m=l m=l SIIITTZ^ 


the other, Q , contains terms linear in A m The situation is different for M odd where 
both solutions take the form of a product of sine functions. See also the discussion in 
[211 for the case when q is a root of unity and the second solution Q 

factorizes into Q~^ and a ’’complete string” polynomial due to a loop algebra symmetry 
of the model |2h| . Here we will remove the associated degeneracies in the spectrum. 

In this article we give an explicit construction of the Q’opsi’a-tors behind the two 
solutions for even M and relate via a limiting procedure the linear terms in the 
second solution to the decomposition in ®. For the construction of Q-operators one in 
general assumes analyticity with respect to the variable z = by requiring that the 
Q-operators commute for arbitrary values of the spectral parameter, {Q{zi),Q{z 2 )\ = 0 
CHI. This usually prevents the appearance of linear terms in A. 

The key to obtain these terms is the introduction of quasi-periodic boundary condi¬ 
tions on the finite lattice taking the limit to periodic boundary conditions at the very 
end of the construction. The twisted boundary conditions depend on a generic twist 
parameter a where the value a = 0 corresponds to ordinary periodic boundary condi¬ 
tions. For a ^ 0 two linearly independent solutions to m without linear terms can 
be constructed and one finds a simplified expression for the six-vertex fusion hierarchy: 
denote by the transfer matrix with spin n/2 in the auxiliary space then one has 


r"“Q+(A + §)Q-(A - I) - g-“Q+(A - f )Q-(A + f ) 

“ q-S^-a _ qa+S^ 


(15) 


Here denotes the total spin operator acting on the chain. By analytic continuation 
of this formula with respect to the spin variable n/2 of the transfer matrix one obtains 
in the limit a —> 0 an alternative expression for 0, 


t(A, C) = lim 


C"Qi{q\hQaiQ^ch - CQtiq\~hQaiQ^&j 

q—S^—a _ qa+S^ 


(16) 
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which reproduces all the desired properties m of the trace function, in particular the 
appearance of linear terms in A, 

q^h{X,C = q^) = \-g{q‘^^)+9{q^^) • (17) 

In Hlh|) we have identified on the right hand side of the equation Qq(A) = Qa{q^) 
and g{x),g{x) in (El) are operator valued polynomials whose degrees in each fixed spin 
sector are given by 


degg\s- = M - |5^| and deggfis;. < M , (18) 


where the upper bound for deg [52 is assumed when ^ 0 but it is strictly smaller 
than M when = 0. This will be derived in Section 4. Before that we present in 
Section 3 our concrete construction of the Q-operators, which will relate the appearance 
of the extra variable C to a special restriction of a Verma module of the upper Borel 
subalgebra Ug{b^) C Uq{sl2)- 

It should be emphasized that our construction of Q-operators for the quasi-periodic 
chain has applications beyond the focus of this article and it would be interesting to see 
whether similar formulae also occur in the treatment of the correlation functions for the 
finite XXZ chain i^Eiini- These approaches rely on the Bethe ansatz and we stress that 
the Bethe ansatz equations are also contained in ED when setting n = 1 and shifting 
A —> A + 1 / 2 , 


M 

[A — Am + 1](J 

m=l 


q-^Q+{\ + l)Q-{\) 

q-S--a 


-q^^QtWQaiX + l) 

_ qo+S^ 


(19) 


In the context of the Liouville model the analogue of this identity has been referred to 
as quantum Wronskian m- Evaluating it at A = — 1 gives the Bethe ansatz 

equations above and below the equator naEu. However, the above equation (dni) is 
of a simpler form than the Bethe ansatz equations and it is desirable to understand its 
consequences in the thermodynamic limit. 

An analogous construction of Q-operators for the eight-vertex or XYZ model is 
currently under investigation m- 


2 The six-vertex fusion hierarchy 


The inhomogeneous six-vertex model on a finite lattice is associated with the quantum 
loop algebra Uq{sl 2 ) in terms of which the transfer matrix and fusion hierarchy can 
be defined. Consider the spin n/2 representation of the hnite quantum subgroup 
Uq{sl2) C Uqish), i.e. 


7rW(e)|A:) = [n - k + l]q[k]q \k - 1) , (20) 

7rW(/)|A:) = |A: + 1), 

7rW(gt)|A:) = q^-’^\k) 
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with /c = 0,1, n. The transfer matrices of the model are built up from the intertwiners 
of the respective evaluation modules. Define 


L = 


z q2+^ — q 2 z{q — q 

(q — q~^)eq~^ z q~i^^ — q^ 


G Uglsh) ® EndD 


( 21 ) 


where V is identified as the two-dimensional representation space of and z = q^^ 
as the spectral parameter. Then the inhomogeneous transfer matrix of spin n/2 is 
introduced by setting^ 

rW(z) = Tr q'^^^^LM(zCl) ■ • • Li(zC?) G EndV^^ . (22) 

ttW 


For the moment the length of the spin-chain M can be any positive integer and a is 
the twist angle. The set {Cm = some unspecified ’’generic” inhomogene¬ 

ity parameters. The above transfer matrices constitute the six-vertex fusion hierarchy 
satisfying the functional equation m 


rW(z9-+i)r«(z) 


M 

ri”+‘>(V) n 

m=l 


M 

l)+4-i)(zg-+2) J]«^-l). (23) 

m=l 


In the homogeneous limit, (m Ij we obtain from the transfer matrix T as logarithmic 
derivative the Hamiltonian for the finite XXZ chain, 


Hxxz 


-(^-O 


Tgjz) 

^ dz ^ {zq^ - l)^ 


M 

--y 
2 ^ 

m=l 




+ O-mO" 


y r^y 


m'-' m+1 




q + q 


-1 




m^m+1 


- 1 ) 


(24) 


Here the twisted boundary conditions manifest itself in the identification 


"M-i-i ~ q '2 


± 

1 • 


The well-known symmetries of the model are expressed in terms of the following com¬ 
mutators 

[rM(z),rW(u;)] = [Ti^Hz),s^] = [rW(z),e] = o, (25) 

and 

mTi^\z) = y^{z)m (26) 

where the respective operators are defined as 


M M M 

m=l m=l m=l 


These symmetries hold for spin-chains of even as well as odd length. For later purposes 
we also compute the value of the transfer matrices at the origin, 




(_)M qi»-snh 


(-) 


Mq 


.n{S^—a) 


q‘ 


S^-a 




(28) 


^Notice that we have changed our conventions in comparison with ESI mi ED and have made the 
following change in notation, T^"\z) —> (aq"). 
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3 The auxiliary matrices 

Our construction follows the main steps of the previous works m and [la 1210011111 ED 
with some minor modifications. We therefore only state the main results and briefly 
comment on the steps involved to accommodate twisted boundary conditions. The 
reader is referred to the aforementioned literature for the technical details. 


3.1 Definition 

For the definition of the auxiliary matrix we employ the following representations vr^ = 
TT^{z;r,s) of the upper Borel subalgebra Uq{b+) C Uq{sl 2 )- Let k G N>o then 


7 r+(eo) |A:) = z\k + l), -k^ {q ^ ) \k) = {q ' 2 -)\k) = rq ^ |A:) , 

1 2k —2k 

7T+{ei)\k) = - - —-—- 1^-1), 7r+(ei) |0) = 0, r,s,z£C (29) 

[q-q 


and set 


TT 


:=7r+ocu with {ei,eo,q ^ ,q ^ ^ {eo,ei,q 2 ^q 2 } 


(30) 


The representation vr is a particular restriction of the representation introduced in EH- 
When q is a primitive root of unity of order N we set N' = N if the order is odd and 
N' = N/2 if it is even. Then the above infinite-dimensional representation becomes 
reducible and can be truncated EH, 


7 r+(eo)|iV'-l) = 0 . 


(31) 


For roots of unity this truncation will always be implicitly understood. In order to 
unburden the formulae we will often drop the explicit dependence on the parameters 
{z,r,s} and the representation 7r+. Set 


( zl:q 2 +^-q 2 {q-q ^eoq 2 ° 

t -i\ -^+1 ^ 

\ [Q ~ Q )^iQ ^ zr q 2 ^ —q 2 


G 17,(5+) 0 End y, 


(32) 


then = (vr'*' 0 1)T is the intertwiner of the tensor product vr'*' 0 Notice that 
the intertwiner for the representation tt~ is obtained via spin reversal, i.e. 


2,^,- = (1 0 <T*)£jr+ (1 0 Cr"^) ■ 


(33) 


Dehne the auxiliary matrix in terms of the intertwiner £ as the trace of the following 
operator product, 

Qa{z]r,s) =TTq°'^^®^iLM{zCM\r,s)-■ ■ SLi{z(\-,r,s) . (34) 

7r+ 

For later purposes we also define the special limits 

Q'^{z) = limQa{0]l,s)~^Qaiz-,l,s) and Q~ = ■ (35) 

s —>0 
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By definition of the auxiliary matrix the operators are well defined and we have for 
roots of unity 


q^ = l: Q„(0;r,s) = 


(36) 


while for 


generic q : Qa(0;r, s) = 


qoi-S^ — 


with < 1 for > 1 . (37) 


After a suitable renormalization (see Section 4.1) the eigenvalues of the operators (13511 
yield the two linearly independent solutions to Baxter’s TQ equation mentioned in the 
introduction. 


3.2 Spin conservation and reversal 

The matrix (1511) commutes by construction with the fusion hierarchy and preserves two 
of the symmetries (15111 [55115^. 

[Qa{z\ r, s), 5^] = [Qaiz; r, s), 6] = 0 . (38) 

Because of spin conservation the dependence of the auxiliary matrix Q on the parameter 
r can be easily extracted, 

Qa{z; r, s) = r^~^"Qa{z- 1, s) = r““‘^"Qa(z; s) (39) 


Without loss of generality we can therefore set r = 1 in order to discuss the behaviour 
of the auxiliary matrix under spin reversal. Employing once more spin conservation we 
easily find from 


(u+|£|u+) 
(u_| £ |u_) 


i+— 

zs q 2 — q 

-I-hi hi 
z q 2 — g 2 


hi 

2 


/ \/ —1 —2 —1 ^—hi hi\ 

= {—zsq){z q s q ^ — q ^ ) 

/ N/ _i _2 —1 —hl\ 

[—zq){z q s sg ^ 2 —q 2 ) 


the identities 


^Qaiz;s,{C^})m.= {-z)’^q’^s 2 ^^Q^{z ^q '^s ^s,{Cj})‘nC 




>2 

sm 


(40) 


and 

iHQo(z, q] s)iH = Q-a{zq^s, q~^]s~^f . (41) 


3.3 Commutation of the auxiliary matrices 

For generic q we can employ the concept of the universal i?-matrix and the fact that 
7 r'’'(2;; r, s) is the restriction of a well-known evaluation Verma module of Uq{sl 2 ) to 
conclude that the intertwiner, say 5, of the tensor product 

7 r“''(2;; r, s) 0 r\ s') 
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exists. Consequently, the Yang-Baxter relation 512 XI 13 XI 23 = '^ 23 '^i 3 'S'i 2 holds and the 
auxiliary matrices commute among each other, 

[Qaiz-,r,s),Qa{w,r',s')] = 0 . (42) 

Here we have implicitly used that 

[£, 7r+(/ii) 0 a^] = [S, Tr~^{hi) 0 7r^(/ii)'] = 0 (43) 

in order to ensure compatibility with the quasi-periodic boundary conditions. The above 
commutation relations are a direct consequence of the intertwining property of £ and 

5. 

When q is a root of unity we cannot use the existence of the universal i?-matrix and 
have to construct the intertwiners S explicitly. This has been done for N = 3,4,6 in 
[13 ED- Numerical checks have been carried out for N = 5,7,8 and the commutation 
relation was found to hold. We therefore assume as a working hypothesis that the 
intertwiner S also exists in the root of unity case for general N. 

From (Bll) it follows that the eigenvalues of the auxiliary matrix are polynomial in 
the parameters z, s and that the eigenvectors of the Q-operator do not depend on them. 

3.4 The TQ equation 

One of the essential properties of the auxiliary matrix Qa is that it satisfies the functional 
equation 

Ta{z)Qa{z;r,s) = 

Qa{zq^;rq~^, sq~‘^) - 1) + Qa{zq~‘^-,rq, sq^) - 1) (44) 

m m 

which one derives from the following non-split exact sequence of representations nznm 

0 — 7r~^{zq‘^;rq~^, sq~‘^) ^ 'k~^{z‘, r, s) 0 {zq~‘^,rq, sq^) —> 0 . (45) 

The proof can be found in EZlEEl, here we have only to incorporate twisted boundary 
conditions. To this end note for instance that the inclusion i \ T:^{zq^',rq~^ ,sq~‘^) ^ 
TT'^{z', r, s) (Si in (??) is given in the following form 

\k) ^ \k + 1) s + Ck\k) S V- 

where is some coefficient whose explicit form is not relevant here. Then one easily 
verifies that 

toTT^{q°‘^^;z,rq-^,sq~'^) = 7r+(g“^i; z, r, s) o z . 

Similarly one shows that the projection p : 7r'''(2:; r, s) 0 —> 7r~^{zq~‘^,rq, sq'^), 

p : \k) (—>■ c'l^ \k) S v+ 

in (??) commutes with the twist operator as well, 

p o Tr~^{q°‘^^; z, r, s) = 7r~^{q°‘^^; z, rq, sq^) op . 
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Setting r = 1 and taking the limit s ^ 0 we find 


Ta{z)Qt{z) = 

m 


l)+q^^^-^^'>Qt{zq-^)\{{zClq^ 

m 


1 ) ■ 


As we will see below, this relation is identical to Baxter’s famous TQ equation mentioned 
in the introduction after an appropriate rescaling of the auxiliary matrix; see Section 
4.1. The eigenvalues of the operators defined in (ESI) then coincide with the two 
linearly independent solutions of m discussed in the introduction. 

As an immediate consequence of the TQ equation and the fusion relation we have 
in terms of eigenvalues the identity 




= (2r")<3t (^9") E 

i=i 


Q+(V€-^)g+(V£-n-2) 


(46) 


This formula is easily proved by induction. 


3.5 Functional equation at roots of nnity 

When g is a root of unity the vital information on the spectrum of the auxiliary matrix 
is encoded in the following functional equation which is a straightforward generalization 
of a previous result m to twisted boundary conditions^, 


Qaizq'^ls] s)Qa{z]t) = 


^S^-a 


Qaizq‘^/s;stq 


-2\ 




■ (47) 


It implies the following decomposition of the eigenvalues 


n± 

Qa{z-,s) = QaiO)Q+{z)Q~{zs), Q^{z) = Z xf) (48) 

i=l 

where the last factors are polynomials with a total of M = n+ + n_ roots and are 
related by the following formula 


N' 


Qa{^)Qa {z) = g(2^'+l)(5" “)g+(^) ^ 


£=1 


Qtizq^QQtizq^^-"^) 


(49) 


Employing the transformation laws EDI), ED) of the auxiliary matrix under spin-reversal 
we deduce that 

M 

n± = yT5". (50) 

Bn order to facilitate the comparison with imrrn note that we have made the following changes. In 
order to match the result in 1221 set s = and r — ■ We have also redefined the fusion hierarchy, 

T^'^\z) —> T^"~^\zq'^), and shifted the parameter r in Idll by r —> rq~^. 
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From the above identity (HU) between the eigenvalues of we now derive the following 
expression for the fusion hierarchy 




(2N'+l+n)iS^-a)Q+f 2n^Q+f ^ ^ -1) 

q Q+(^^2.)Q+(,g2.-2) 


q 


{2N'+l+n){S^-a)^+ 


Ql{z)Qt{zq 


2n\ 


N' 


J-n+l J 


After inserting the value for (5o(0) we obtain 


i qS^-a_qa-S^ 


(51) 


As long as a 7 ^ OmodA^'jS'^ this expression holds true for M,N even and odd. In 
particular we have for n = 1 the identity 


M 


TTn ^^2 _q^’' '^Qt{zq)Qa{zq ^) - q"^ ^’'Qtizq ^)Qa{zq) 

XX U ^SmQ) 


rn=l 


qS^-a _ qOL-S- 


which has been called ’’Quantum Wronskian” in the literature m- 


3.6 The algebraic Bethe ansatz for generic q 

In [20] the spectrum of auxiliary matrices constructed in jOZj have been investigated 
on the basis of the algebraic Bethe ansatz. The results in m apply to the inhomoge¬ 
neous XXZ-chain with twisted boundary conditions. Denote the monodromy matrices 
associated with the transfer and auxiliary matrix by 

T=(7rW®l)(?“'^®iLM---Ti= ( 


and 

Q = (vr+ ® = {Qki)k,e>o, Qki = {k\ Q \i), 


respectively. If the quantum space Ti carries a highest (or lowest) weight representation 
of the quantum group with highest weight vector |0).^ then the eigenvectors and eigen¬ 
values of the Q-operator can be computed via the algebraic Bethe ansatz. Setting r = 1 
in (P|) and denoting by = {l,e'\2,\k,e) the matrix elements of the T-operator one 
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finds 


n+ 

Qaiz)'[lB{xp | 0 )„ = 

i=i 

g (0|Qfcfc(^)|0)^ n (1 _ zxpq-^^^){l - zxpq-^^-^) 


n+ 


ne(.<)|o> 


H 


i=i 

n+ 


i=i 


where are the Bethe roots above the equator. This formula has been proved directly 
for n+ = 1,2,3 in |3()j and further checked for consistency for arbitrary n+. In the 
present case of the XXZ chain we have 

|0)„ = u+® n = {C^)®^ (53) 

and 

mkkizMn = - 1 ) • ( 54 ) 

m 

The above expression from the algebraic Bethe ansatz then yields for the case of generic 
q the same information on the spectrum which we had previously derived from the 
functional equation at a root of unity, 


n+ 


Qa{z)'[lB{xj) | 0 )„ 

i=i 


^S^-a 


Qtiz) Qt{zs)Y^ 




n+ 


k>0 


Qt izsq-^’‘)Qt {zsq-^^-^) 


i=i 


where we have set as before 


n+ 


A// 

Qtiz) = - zxj), n+ = — -S^ 

j=i 


(55) 


Thus, the analogous result concerning the decomposition of the eigenvalues into 
applies also here. Notice that for generic q one has to put further restrictions on the 
twist parameter a m, 

|^-a±M/2| ^ |^|±1 > 


in order to ensure absolute convergence. Analytically continuing from this region the 
relation m holds also true for generic q. 

An alternative to the Bethe ansatz is to follow the analogous line of argument as 
presented in m- Note that when we set s = in (EHl) it follows that 

s = : ei |n) = 0 . (56) 
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Hence, the infinite-dimensional module tt'*' splits into a finite n-dimensional part VH<n 
given by the linear span of the basis vectors {|^)}fcZQ and an infinite-dimensional space 
W>ri which is the linear span of the vectors {|A:)}^^. Note that W>ri is left invariant 
under the action of Uq{b^). Restricting the Xl-operator onto these two spaces, one finds 

{7r+(2; 1, s = ® l} £|vk<„ = (1 ® ^® l} L{zq^){l 0 q^^^) (57) 

and 

{7r+(z; 1, s = g2-) ® i} = (i ® {n+izq^^; 1, s = 0 l} £(1 0 . 

(58) 

Employing these identities together with the invariance of W>n one can split the trace 
in the definition of the auxiliary matrix into two parts leading to 

Qo^iz- S = g2-) = ^ ^ ^-2n) _ ( 59 ) 

Apart from the factorization of Qa into this corresponds to the identity (|^ . As 
it turns out this result is already sufficient for our present purposes, i.e. the analytic 
continuation of the fusion hierarchy in the spin-variable. See equation m below. 


4 The trace function 

In order to make contact with the trace function used in m we now reparametrize the 
fusion hierarchy and the Q-operator and take M to be even. 


4.1 Reparametrization 

Henceforth we set z = q'^^, Cm = 9”^"’ 


and dehne the rescaled fusion hierarchy as 

^ {zq)-fTi''\z) 1 


(60) 


This renormalization corresponds to the following choice of the six-vertex R-matrix 
which is in accordance with the conventions used in na. 


r(A) 


/ [A + l]g 0 0 0 \ 

0 [A], 1 0 

0 1 [A]g 0 

\ 0 0 0 [A + l]g / 


(61) 


Thus, the fusion hierarchy is now expressed as 

tM(A + !!±i)(„(A) = (l;‘+‘l(A + I) p [A - A„ + 1], + (l,”-‘)(A + eM) p[|A - A„], (62) 


m 


m 


with the quantum determinant being 


t!.»>(A)=t<“>(A) = n[A-A„ + il, . (63) 

m 
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Let us now turn to the re-definition of the auxiliary matrix. With respect to the decom¬ 
position 

£t 


£(A) = 




the matrix entries are now chosen as 


*9 / \ I '‘i+i 

a+ _ A-i C A 2 - g ^ 2 

q-q ^ 

I 1 \ ^Q + l 

£+ = ( 2 q ^eoq 2 , 

>11-1 

= C ^ Q Q 2 ei, 

_ ._i r q 2 — q 2 

L_ = C ^ — 


s = e,z = q 


,2A 


q-q 


-1 


Notice that this re-definition corresponds to the overall scaling factor 

Qaiz] r,s = C^) ^ Q„(A; r, C) := , ^ 0 

\Q Q ) m 


(64) 


(65) 


and the eigenvalues of the auxiliary matrix in each spin-sector have now the following 
decomposition 

Wx_m ^a-S^X'S^n+l 


Qa{X; r,C = qn=^a q^ ^ Qt{X)Q-{X + A') 


( 66 ) 


with 


and 


= 


q 


-N’(S^+a) i^'(a+S^)]g ■fN_i 
[a+s-], ’ nq - L 

^ if q generic 


n± 




± -2£^ 
xf = q 


i=l 


Here we have used the sum rule 


1 lla^i lll'i IKm g„{0,<l) -,a+S* 


L L 

2=1 2=1 m 


(67) 


which follows from combining dm) with m and setting z = 0. Employing this decom¬ 
position the TQ equation is equivalent to 


M 


M 


ta{X)QtiX) = q-^^QtiX + 1) n + ^"QtiX - 1) n + !]<? • (68) 


m=l m=l 

The expression for the transfer matrices of the fusion hierarchy is now 


-1),,, _ <i-"‘Qi(x + i)Q^(A - j) - q"‘QH\ - ;)g^(A + 


e-‘’(A) = 


q—S^—a _ qu+s^ 


, n G N . (69) 


It is the last expression respectively m which we want to analytically continue in n in 
order to obtain the trace-function used in the description of the correlation functions of 
the infinite XXZ-chain. 
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4.2 Analytic continuation and the limit a —0 

As mentioned in the introduction we now analytically continue the expression for the 
fusion hierarchy (EU, (isni) in the spin-variable n/2 by defining the following operator 


T^{zX) 


(-) 


qS^—a _ qCx—S^ 


C“ ^'S = C^)-C'^" "Qa«;S = C 


qS^—OL _ qOL—S^ 


Qa(0) 


-1 


(70) 

(71) 


Note that this analytic continuation is unambiguous as the operators on the right hand 
side of dZU have only a polynomial dependence on the spectral parameter z by construc¬ 
tion. The rescaled counterpart of dZU gives the result m stated in the introduction, 


*a(A,C) 


Ta{z = q^\C) 

qM{X+l/2)^q_q-l^M H 


(72) 


First note that we recover the fusion hierarchy m respectively dsni) when setting Q = 


T^{zX = q^)=T^^-^\z) . 


(73) 


Now, as an easy example for the occurrence of terms linear in A in the matrix elements 
we evaluate dZOl) at the origin z = 0. Then we have for nonzero spin ^ 0 


^-a+S^ _ ^a-S^ 

(-)^ limTa(0, C) = lim ^^ — - 

^ ' a^o ^ ^ a^o q-°‘+S - q°^-S^ 


q 




q 


-s^ 


(74) 


For vanishing spin = 0 we set C = 9'^ obtain 


(-)^ lim r,(0,C = <7^) = A. 


(75) 


These last two relations correspond to the defining equations (5.3) in jTl], see also m in 
the introduction of this article. The occurrence of linear terms A in the matrix elements 
is not restricted to the zero spin sector but occurs more generally. To see this we now 
derive the analogue of Lemma 5.1 in [El- 

Lemma. For M even the analytically continued fusion hierarchy decomposes 
in the limit a —> 0 into a sum 


lim (z, C = g^) = A • G(z) + G{z) (76) 

a^O 

where the operators G, G are polynomial in the spectral variable z and in each fixed spin 
sector ^ Q have degrees 

degGIs^ = M — and deg G\s^ = M . (77) 

If = 0 then 

degG| 52 =o = M and degGl^z < M . (78) 

According to the rescaling this obviously implies 03). 
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Remark 1. For comparison with Lemma 5.1 in m we have to identify z ^ Ci aiid 
M = 2p — 2, compare with 0 in this article. Moreover, we do consider here the degree 
of an operator in a whole spin-sector not a single matrix element as in mi. 

Remark 2. Notice that for M odd the linear terms in A are absent, since the terms 
containing Tr^+ which contains a pole in the limit a —> 0 can never occur in a 
matrix element of the Q-operator. See the proof below for an explanation. 

Note further, at roots of unity and periodic boundary conditions a = 0 the linear 
terms have not been observed in m and This is explained by the fact that 

the root of unity limit does not commute with the limit a — 0 . As long as a 7 ^ 0 
(and generic) the root of unity symmetries discussed in (21112311211221 are not present, 
whence we arrive here at a different result. 


Proof. Since we are only interested in determining the maximal degree of the 
respective polynomials in z we can set Cm = 1 without loss of generality. As an additional 
preparatory step we need to make contact with the Casimir operator which has also been 
used to define the trace function in mi. Set r = 1 in (| 2 ^. then we have 


[ei,eo] = z 


— g 

g - g-i 


(79) 


and 


gqhi 1 _|_ g ^1+1 1 _)_ 


(80) 


This identity corresponds to equation (5.4) in mi) s®® m in this article. Now consider 
a general matrix element of the auxiliary matrix 


g«(z;5)f' = Trg“"i®iT: 


£m 


0^1 
' .^£ 1 ) 


Em — Z] — 25"^ 


(81) 


Since the Q-operator preserves the total spin the matrix elements always occur in 
pairs and we compute 


Tr{X = zg Tr{A (g - g"^)^^: - 

7r+ 7r+ 

= zg TV{A (l + s-g-'^i^^-sg^i^^)} (82) 

7r+ \ / 


The numbers of the various operators occurring in a matrix element is given by 


+ = = n^ and #£+ = #£;= ^ 

" m 


1 






(83) 


Using the above relations any matrix element of the Q-operator in a fixed spin sector 
can be expressed as a linear combination of terms of the form 


g„(zC; s = c-')|' = {zqT^^^' (c"' - 9”^) 


h. \ n--n- -/ 


zq 


1 _ ii hi\ 


2 — (7 2 


n c^Q 

i=l 
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where the coefficients do not depend on z and Oi = 0 , ±1 can vary in each factor, 
but there can be terms present for which all the Uj’s are equal. To see this note that 
= 0, while the commutation of with only produces powers in q. The 
operators and commute according to (jT^. Let us distinguish the two cases 

Choose a matrix element with = 0, i.e. Em = £m- This is possible in all 
spin-sectors. Then the term of maximal degree has the coefficient 


Q„(zC;s = C ^)| = 


Tr 9 


‘(c-‘ 

x(c 


1 + -in \ 2 
zq ^ 2 — g 2 


2 — q ^ 2 


-|5- 


. 




zq^ 2 — g 2 

7r+ 

Here we have set a = sgn . From this we infer using m that 


°‘Qa{zC;s = C '^)-C ^‘QaizC '-;s = C) 


—2 


a-S^ 


— 1 . 


(gS-a_^a-S^)Q^(0) 


(zq) 


M 


—S^-a 


-c 


a+S^ 


(q- 


S^-c 


r,a-S^ 


)Q«(0) 


ga{l+S^)hi _|_ 


Setting (^ = q^ we conclude similarly to our previous calculation at z = 0 that the linear 
terms in A originate from coefficients containing Tr^+ which develops a pole in the 
limit a ^ 0. Thus, we deduce that degG = M and degG < M when = 0. When 

7 ^ 0, on the other hand, we have as degrees degG = M and degG = M — IS"^!. Note 
that for this it is crucial that 2|S'^| is an even integer which is only the case when M is 
even. 

Now choose a matrix element with 7 ^ 0. This is always possible as long as 
|5^| < M/2. All we need to show is that the just derived degrees are not exceeded. The 
term of maximal degree has now coefficients of the form 


Qa{zC,\s = C ^)i = {zq)^ Trg 


(a+S^)hi 


n 

2 = 1 


7ig-0-ihl _|_ 


As in the earlier examples we need to determine the maximal degree of the term which 
contains Tr.„-+ yielding the linear dependence on A in the limit a —> 0. The degree 
of G is maximized by choosing a matrix element with n^g/ = \S^\, where the coefficient 
of the term with all cTj = sgn is non-vanishing. Then we have 


r -c 


a-S= 


= const. - C)Trq<^’^^ + 


from which we infer in the limit a —> 0 that the maximal degree of G in a fixed spin 
sector is again M — |5^|. For the remainder polynomial G we hnd as before that its 
degree is strictly smaller than M if = 0. □ 
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4.2.1 Example: M = 4 and = 1 

Let us consider a simple example for the homogeneous chain to show how the linear 
term in A emerges in a matrix element of the operator (1701) . Setting M = 4 and = 1 
we choose the matrix element 

Qa{z-,s)-tXX = -q-^fzq(l + s-q-^^-^ 

7 r+ V / 

= zqTtq‘^^^{z^s^q^+'^^ - 2zsq + q-’^^) (l + s - - sq^^+^) 

= -zsq‘^{z‘^s + l + 2z{l + s)}Trq'^’^^ + ... (84) 


In the last line we have only written out the terms which will give rise to the linear 
dependence in A. Namely, inserting this expression into dZOl) we find 


Ta{z-,C) 


-+++ 

+-++ 


C-^Qo.izC;C')-C-^QaizC-^;n 

Qa(0)(gi-“ 

-zq'^iz'^ + l + 2z{C + C^)) 


a _ 


-zq\z^ + 1) 


a _ 




Trg"^i + ... 




Trg“'^i + ... 


(85) 


Let us distinguish the case when q is generic and when it is a root of unity. If q^ 
then 


Qa{0) 


1-q 


-2N'a 


^OL—1 


7 I—a 


and 


7r+ 


l-q 


-2N'a 


qoi _ q-a 


= 1 

( 86 ) 


On the other hand we have for generic q 


Qa{0) 


q 


a—1 



and 



1 

q°‘ - q-° ’ 


(87) 


where these expressions are understood as analytic continuation from the region where 
the trace converges. Thus, setting (" = we arrive at 


lim T„( 2 ;;C) 

a —>^0 


+-++ 


-zq'^iz'^+ l + 2z(q^+ q~^)) lim ^- - -h ... 

^ VI _ q-a 

-X zq^iz'^+ l + 2z{q^+ q-^))+ ... (88) 


where the coefficient of A is of degree M — | | = 3 in 2 in accordance with our lemma. 


5 Conclusions 

In this paper we continued a previous study EHl nil 1301121 ED on the explicit construc¬ 
tion of operator-solutions to Baxter’s TQ equation m- In terms of eigenvalues this 
equation is a second order difference equation and its theory resembles closely the one 
of second order ordinary differential equations. In the case of the XXZ chain we dis¬ 
cussed the existence of two linearly independent solutions when quasi-periodic boundary 
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conditions are imposed by explicitly constructing the relevant Q-operators using repre¬ 
sentation theory. Because it is usually required that Q-operators should commute for 
arbitrary values of the spectral parameter, [Q{z)^Q{wy\ = 0, this precludes by choice 
of the construction method the possibility of obtaining “non-analytic” solutions to the 
TQ equation: that is, solutions which obtain the linear terms discussed in the text and 
which have been postulated in m for even chains with periodic boundary conditions 
and when q is generic. The result of this paper is that such solutions can arise by 
taking the limit from quasi-periodic to periodic boundary conditions in the explicitly 
constructed Q-operators. Recall that for generic q this limit required a careful analysis. 
First one had to choose the twist parameter a such that convergence of the trace over 
the infinite-dimensional auxiliary space is guaranteed. In a second step we then analyti¬ 
cally continued the matrix elements in a from the region of convergence to the complex 
plane which enabled us in the final step to discuss the limit a —> 0. To complete the 
investigation by computing the spectra of the Q-operators in this limit one would need 
to know the explicit dependence of the Bethe roots on the twist parameter. This is a 
rather formidable challenge as the solutions to the Bethe ansatz equations are in general 
not known. 

The main motivation for our construction of this “non-analytic” Q-operator has been 
the relation with the recent developments in the computation of correlation functions as 
explained in detail in the introduction. The alternative expression m for the special 
trace of the monodromy matrix © entering the ansatz in m shows that there is a 
more fundamental quantity in which the correlation functions can be expressed and 
provides a different point of view on the role of the trace function In future work 
it needs to be explored whether there are concrete practical implications of the identity 
dEI) which facilitate the computation of correlation functions. For instance, one might 
ask whether one can insert each of the two terms in the difference dm separately in 
the ansatz for the correlation functions and if they satisfy identities analogous to the 
quantum Knizhnik-Zamolodchikov equations. Of course one has to keep in mind that in 
order to perform the limit a ^ 0 both terms will be needed at the end. However, such 
an investigation might yield further insight into the analytic structure. 

Note added in Proof: After this paper was completed the work m appeared 
which extends the investigation of correlation functions to the inhomogeneous eight- 
vertex or XYZ model. An analogous Q-operator for this model would be helpful as it 
would simplify the computation of the analogue of the trace function over the Sklyanin 
algebra. The present discussion provides a first step towards this aim. 
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